Global spin alignment of vector mesons and strong force fields in heavy-ion collisions 
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Introduction. Particles and fields are two fundamen- 
tal forms of matter in our natural world. Particles are 
excitations or quanta of their corresponding fields, which 
are characterized by quantum numbers such as mass, 
charge, spin, and parity. Photons, for example, are 
quanta of electromagnetic fields that mediate electromag- 
netic forces among charged particles. In collisions at high 
energy scales, gluons appear as quanta of color fields that 
mediate interaction among quarks as elementary parti- 
cles of strong interaction. At low energy scales, strong 
interaction is often characterized by mesons as effective 
degrees of freedom of quark and gluons, whose existence 
was proposed by Yukawa in 1935 in analogy with elec- 
tromagnetic fields [I]. We now know nuclear forces as 
strong interaction at low energies have many components 
of meson fields [2] B] c, 7, p, w, etc.. As the energy scale 
increases in nuclear reactions other meson fields carrying 
strangeness quantum number may come into play, such 
as K, ¢, etc.. As quanta of strong force fields, all these 
mesons with their specified quantum numbers were dis- 
covered before 1960s. Due to their short Compton wave- 
lengths, experimentally it is much easier to detect par- 
ticles than wave-like fields that are elusive and strongly 
fluctuate. 


Recently, for the first time, the print of one kind of 
strong force field could be detected through the global 
spin alignment of vector mesons in heavy-ion collisions in 
the STAR experiment [4]. This is another breakthrough 
after the one in the measurement of the global spin po- 
larization of A in heavy-ion collisions [5]. The global 
spin polarization arises from the initial orbital angular 
momentum (OAM) in nuclear collisions that is partially 
converted to the local OAM or vorticity leading to the po- 
larization of hadrons through spin-orbit coupling in the 
interaction [6}8]. The spin alignment of vector mesons 
in heavy-ion collisions along the OAM direction was first 
proposed by Liang and Wang some years ago [7]. The ob- 
servable is the 00 element poo of the spin density matrix 
for vector mesons which can be measured by the angu- 
lar distribution of their strong decay daughters in their 
rest frame [9]. The STAR data for po9 show a surpris- 
ingly large positive deviation from 1/3, which is orders of 
magnitude larger than predictions by conventional mech- 
anisms [4]. In Ref. [10], Sheng, Oliva and Wang proposed 
that the local correlation or fluctuation of the @ meson 
field can produce a large positive deviation for poo from 


1/3 and thus provided a possible explanation. Such ¢ 
field may originate from non-perturbative strong interac- 
tion coupled to s and 5 and is connected with vacuum 
properties of quantum chromodynamics (QCD) [I] [12]. 
Vector meson spin alignment The spin state of a sys- 
tem of spin-S particles can be described by the spin 
density matrix ô, which is a (25 + 1) x (28 +1) com- 
plex and Hermitian matrix for spin-S particles with unit 
trace. The number of independent real variables in 6 
is 45(S +1). For example, for spin-1/2 particles, there 
are 3 independent real variables corresponding to a po- 
larization vector P as in p = (1/2)(1 + ø - P) where 
o are Pauli matrices. For spin-1 particles such as vec- 
tor mesons, there are 8 independent real variables cor- 
responding to a polarization vector P (3 variables) and 
a symmetric traceless tensor T;; (also called the tensor 
polarization, 5 variables). As a spin-1/2 particle, the A 
hyperon’s spin polarization can be measured through its 
parity-violating weak decay A — pz, since the preferred 
direction of the daughter proton’s momentum is along its 
spin in the rest frame. However, this is not the case for 
vector mesons since they mainly decay through strong 
interaction in which the parity is conserved. So the ele- 
ments of the spin density matrix that can be measured 
are its tensor components T;j and poo is associated with 
them. For the ¢ meson’s strong decay ¢ > K*K~—, the 
daughter particle’s polar angular distribution is given by 


dN 3 
deoso 4 [1 — poo + (3p00 — 1)cos*8] . (1) 


We see that poo appears in the coefficient of cos? which 
can be determined by measuring the polar angle distribu- 
tion. If poo = 1/3, the distribution is a constant indicat- 
ing that the daughter particles are emitted isotropically. 
If poo # 1/3, the probabilities of the vector meson in three 
spin states are not equal and then the emission of daugh- 
ter particles is anisotropic. For poo 2 1/3, the polar angle 
distribution is in a cigar/discus shape (as shown in Fig. 
E. Correspondingly, if poo > 1/3, the vector meson has 
a larger probability to be in the spin state À = 0 so that 
its average polarization vector (not the spin) is aligned 
to the spin quantization axis. If pag < 1/3, the vector 
meson has a larger probability to be in the spin states 
A = +1 and its average polarization vector is aligned 
in the plane perpendicular to the spin quantization axis. 
The quantity poo is thus referred to the spin alignment 
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Figure 1: (Color online) The shape of dN /dcos@ for the daugh- 
ter particle in the vector meson’s decay with poo = 1/2 > 1/3 
(left panel) and poo = 1/6 < 1/3 (right panel). 


of the vector meson. 

Experimental results. In Ref. [4], ¢ and K*° mesons 
are observed by paring of their decay daughters K= and 
Kr, respectively, with subtraction of the combinatorial 
background and application of the so-called event mix- 
ing and rotation methods. Detail studies show that both 
techniques can effectively break the correlation between 
pairs in real events, and the yields of ¢ and K*° from 
two techniques are consistent [4]. The spin quantiza- 
tion direction is chosen to be the normal direction of the 
second-order event plane constructed from charge par- 
ticles collected by the STAR Time Projection Chamber 
(TPC). Then the polar angle distribution of Eq. (1) is an- 
alyzed, and poo is extracted after correction for detection 
efficiency and acceptance. The spin quantization direc- 
tion can be constructed by different detectors such as the 
shower maximum detector and the beam-beam counter, 
and they all give consistent results on the global spin 
alignment of ¢ and K*° M]. In the following, we only 
discuss results with respect to the TPC event plane. 

Figure [2] shows the poo data for ¢ and K*° mesons in 
Au+Au collisions at \/syn— 11.5 to 200 GeV measured 
by the STAR collaboration [4] from a dedicated Beam 
Energy Scan program and multiple years of high statis- 
tics event collection of Au+Au collisions at 200 GeV. The 
centrality categorizes events according to the observed 
number of tracks in each collision. Here the 0% centrality 
corresponds to exactly head-on collisions, which produce 
the most tracks, while the 100% centrality corresponds 
to barely glancing collisions, which produce the fewest 
tracks. The STAR measurements presented in Fig. are 
for the centrality interval of 20% to 60%, where one ex- 
pects the largest OAM among collisions and better signal 
to noise ratio in the experimental analysis. A complete 
set of results for centrality and transverse momentum 
dependence can be found in Ref. [4]. The STAR’s data 


give 


poo = 0.3512 + 0.0017(stat.) + 0.0017(syst.), 
ok” = 0.3356 + 0.0034(stat.) + 0.0043(syst.), 


which are obtained by averaging over results at energies 
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Figure 2: (Color online) Measurements of poo with respect 


to the second-order event plane in high energy heavy-ion col- 
lisions. Stars represent the data for ¢ mesons [4] [13], circles 
represent the data for K*° mesons [4] [I3]. The solid red line 
is a fit to data by the model based on strong force fields [IO]. 


from 11.5 to 62.4 GeV for ¢, and from 11.5 to 54.4 GeV 
for K*° [Ø]. Taking the total uncertainties as the sum 
in quadrature of statistical and systematic uncertainties, 
the results suggest that Bos is above 1/3 with a signifi- 
cance of 7.4 a, indicating a significant global spin align- 
ment for the @ meson. The values of pes however, are 
consistent with 1/3. Measurements in Pb+Pb collisions 
at \/SnN= 2.76 TeV [13] in the pr bin nearest to the mean 
pr for STAR data in the range from 1.0 to 5.0 GeV/c are 
also shown for comparison. They are consistent with the 
STAR data where magnitude of the spin alignment for 
both ¢ and K* lie at 1/3 at the highest beam energies 
within large uncertainties. 

Quark coalescence model. The quark coalescence 
model (QCM) is the convenient theoretical tool to de- 
scribe the hadron production or hadronization in heavy- 
ion collisions. Therefore one can calculate the spin den- 
sity matrix for vector mesons in terms of those of the 
quarks and antiquarks in QCM [A00]. A non-relativistic 
QCM is the most simple and transparent one to this 
end [10], [14]. In QCM, the meson’s spin density oper- 
ator can be constructed as Ôm = pq © fg, where fg and 
fq are spin density operators of the quark and antiquark 
in spin and momentum space respectively. The elements 
Nr. (x,p) can be obtained by projecting pm onto two 
mesons’ spin and momentum states and then taking a 
Fourier transformation with respect to the relative mo- 
mentum of two states. The poo for the ¢ meson in phase 
space reads [10] 


2 72 
ies p) ~ 37 3 (PY (x1, p1) P? (x2, P2)) 


where x; = x + Ax/2, x9 = x — Ax/2, pı = p/2 + Ap, 
P2 = p/2— Ap, the average is taken over Ax and Ap 
weighted by the ¢ meson’s wave function, P(x), p1) and 
Pz(x2, p2) are spin polarization vectors of the quark and 
antiquark respectively, and the spin quantization direc- 
tion is taken along the y direction. 


Equation clearly shows that the vector meson’s 
spin alignment is determined by the local correlation of 
quark’s and antiquark’s polarization functions P,(x1, p1) 
and Pz(x2, p2) inside the phase space limited by the me- 
son’s wave function. 


Local correlation of @ fields. It is well known that 
fermions such as quarks at rest have magnetic moments 
proportional to their spins that are polarized along the 
direction of the magnetic field. For moving fermions, 
the electric field comes into play in the form of p x E 
(the spin-orbit coupling or spin-Hall effect). Similarly 
vorticity fields and vector meson fields can also polarize 
quarks and antiquarks. For s and 5 that make the ¢ 
meson, the spin polarization vector is 


1 
Paa = 3” + 5 exp 
Io 96 
2E a E 3 
nt ne 


where T is the effective temperature of the quark matter 
when s and 5 combine into the ġ meson, € and w are the 
electric and magnetic part of the thermal vorticity tensor, 
and Eg and Bg are the electric and magnetic part of the 
@ field, respectively. Here we have neglected the effects 
from electromagnetic fields since they die away quickly 
in the late stage of the matter evolution in heavy-ion 
collisions. The effects from vorticity fields can also be 
neglected because their magnitudes from the measured 
A polarization are too small to account for the observed 
Pons but their terms are kept in Eq. just as a con- 
trast to those of @ fields. The main difference between 
terms of vorticity fields and ¢ fields in Eq. (3) is the sign 
for antiquarks: vorticity fields do not distinguish quarks 
from antiquarks while ¢ fields do (same as electromag- 
netic fields). We see in the following that this is essential 
to obtained a large vector meson spin alignment poo in 
the coalescence of s and 5 into the ¢ meson. 


Substituting Eq. into Eq. (2), one obtains that péo 
depends on the local correlation of ¢ fields inside the ¢ 
meson’s wave function. A more rigorous approach based 
on the relativistic QCM has been formulated using the 
closed-time-path (Schwinger-Keldysh) method in quan- 
tum field theory [I5] [16]. The result obtained for p% 


where the fields with primes are in the ¢@ meson’s rest 
frame, and Cı and C2 are functions of ms (strange 
quark mass) and mg ( meson mass). We see that all 
terms appear in fields squared. The momentum de- 
pendence of pes can be obtained by rewriting Peo in 
terms of fields in the lab frame using Lorentz trans- 
formation. By taking averages over space-time on the 
hadronization hyper-surface of the @ meson, one obtains 
pen as functions of momentum that can be compared 
with STAR’s data [4]. The parameters are in the form 
of (o3E3/T?) and (93B3/T?) and reflect local fluctu- 
ations of ¢ fields [15]. If we assume that local field fluc- 
tuations are different in the transverse (labeled as 7) and 
longitudinal direction (labeled as z) with respect to the 
beam direction z in heavy-ion collisions, then peo depends 
on two parameters F} = (g4 B3,/T°) = (g3E3 ,/T*) and 
F? =o B3, ft) = (g3 Eg, / T°) The values of two pa- 
rameters are determined by fitting STAR’s data on poo 
as functions of collision energies in Fig. With fitted 
values of two parameters one can predict the transverse 
momentum spectra of a at all available collision ener- 
gies which are in good agreement with STAR’s data [4]. 

Summary and outlook An unexpected large global 
spin alignment of ¢ mesons has been observed by the 
STAR Collaboration [4] in relativistic heavy-ion colli- 
sions. By using the quark coalescence model for hadron 
production, Refs. [I0] [15] [16] provide a good interpreta- 
tion of the data [4]. According to this interpretation, such 
a large global spin alignment of ¢ mesons is induced by 
the local correlation or fluctuation of strong force fields 
on the hadronization hyper-surface. The average values 
(E%) and (B%) reflect the local fluctuation of ¢ fields 
and are expected to be calculable using lattice QCD. 
These studies open a new avenue to investigate prop- 
erties of the strongly interacting quark matter as well as 
non-perturbative properties of strong interaction. In ad- 
dition, relativistic heavy-ion collisions are usually called 
“small bang” in contrast to the Big Bang of the universe. 
In such an analog, hadronization on the freeze-out hyper- 
surface corresponds to the stage of the early universe in 
which particles are decoupled from the interaction dur- 
ing the Big Bang. The vector meson’s spin alignment is 
similar to polarization modes of cosmic microwave back- 
ground radiation. 

We note that such an explanation is still subject to 


debate and further verification. More systematical stud- 
ies in both experiments and theories need to be done to 
clarify the deep physics behind the phenomena. 
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